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ABSTRACT 

We use multi-band photometry to refine estimates for the planetary radius 
and orbital inclination of the transiting planet system HD 209458. We gathered 
1066 spectra over four distinct transits with the STIS spectrometer on the Hubble 
Space Telescope using two gratings with a resolution R = 1500 and a combined 
wavelength range of 290 — 1030 nm. We divide the spectra into ten spectropho- 
tometric bandpasses, five for each grating, of equal wavelength span within each 
grating, and fit a transit curve over all bandpasses simultaneously. In our fit 
we use theoretical values for the stellar limb-darkening to further constrain the 
planetary radius. We find that the radius of HD 209458b is 1.320 ± 0.025 i?j up , 
which is a factor of two more precise than current estimates. We also obtain 
improved estimates for the orbital period P and time of center of transit Tc- 
Although in principle the photon-limited precision of the STIS data should allow 
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us to measure the timing of individual transits to a precision of 2 — 7 s, we find 
that uncertainties in the stellar limb-darkening coefficients and residual noise in 
the data degrade these measurements to a typical precision of ±14 s. Within this 
level of error, we find no significant variations in the timing of the eight events 
examined in this work. 

Subject headings: binaries: eclipsing — planetary systems — stars: individual 
(HD 209458) — techniques: photometric 



Introduction 



HD 209458b is the best-studied transiting planet to date, due in large part to its 
proximity and the resultant apparent brightness (V = 7.6) of its parent star. Although 
HD 209458b's physical characterist ics have been measured more accurately than those of 
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the other eight transiting planets (ICharbonneau et al.l l2000l: iHenrv et al.l 12000 ; |Jha et al. 



2005 



Winn et al 



20051 ). it is the only transiting 



planet whose radius is not consistent with the predicted values for irradiated hot Jupiters. 
Guillot et al.l (119961 ) predicted that hot Jupiters like HD 209458b, which orbit at typical 
separations of 0.05 AU from their parent stars, would have radii that are inflated relative to 
Jupiter. The radius of HD 209458b w as 20% larger than predicted, however, even includ- 



ing the effects of increased irradiation (jBodenheimer et al 
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Bodenheimer et al.ll2003l ; iBaraffe et al.ll2003l ; iLaughlin et al.ll2005al ). In contrast, the radii of 
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The precise determination of the planet's radius is also of critical importance to cor 
rectly interpret the results of several follow-up observations of the planet ary atmosphere (e.g. 



Charbonneau et al 



2002 



Deming et al 



2006af), notably through transmission spect roscopy (Charbonneau et al. 



2005aj) 



upper limits on reflected light (IRowe et al.l 120061) . and the 
study of the planetary thermal e mission through secondary eclipse monitoring (iDeming et al. 
2005bl : iRichardson et al.l l2003al lbl) . 



1.1. Proposed Models for HD 209458b 



Several models have been proposed to explain the s ource of additional in t ernal energy 
required to reproduce the observed size of HD 209458b. I Showman fc Guillotl (120021 ) argue 
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that the intense radiation from the star might create strong winds in the planet's atmo- 
sphere, which would transport energy downward and heat the planet's interior to higher 
temperatures than would otherwise be expected. They calculate that 1% of the incident 
stellar radiation would have to be converted into winds to explain the observed radius of 
HD 209458b. The majority of the transiting planets that have been found to date do not have 
radii that are significantly larger than the predicted values for irradiated hot Jupiters, how- 
ever, which indicates that either (1) winds of this strength do not occur on all hot Jupiters 
or (2) the size of the cores in the other (normal) planets must be inc reased to keep thei r 
radii the same even with the additional energy input from these winds ( Guillot et al.ll2006l ). 
If the first explanation is correct, there must be some kind of distinguishing characteristic 
that would allow significant winds on some hot Jupiters but not others. If the second expla- 

d need significantly larger solid cores than currently 



nation is correct, most hot Jupiters wou 
predicted by current formation models (IGuillot et al.l 120061 ). 



Bodenheimer et al.l (120011 . 120031 ) propose that the req uired energy may be provi ded by 



dissipation from tidal circular ization of an eccentric orbit. iBodenheimer et al.l (120031 ) calcu- 
late that an average eccentricity of 0.03 would explain HD 209458b's enlarged radius if the 
tidal quality factor Q « 2.5 x 10 5 . This is co mparable to the tidal quality factor for Jupiter , 
which is estima t ed to b e between 10 5 and 10 6 ( IGoldreich k, Soterlll966l ; lYoder &: Pealdll98ll ). 
Laughlin et al.l (l2005bl ) measure an eccentricity of 0.018 ± 0.009 using radial velocity data 
from 2004/ 2005, and use Monte C arlo simulations to show that this value is in fact consistent 
with zero. iDeming et al.l (l2005bl ) also found no evidence for a non-zero eccentricity when 
they observed the secondary eclipse with the Spitzer Space Telescope, constraining e sin u to 
be less than 0.0015, where e is the eccentricity and u is the longitude of periastron. Thus 
the planet cannot have a significantly non-zero eccentricity unless the semi-major axis is 
precisely aligned with our line of sight, and even then it is still significantly constrained by 
the radial velocity data. 

Although these results indicate that the planet is unlikely to have a significant constant 
non-zero eccentricity, it is possible that the eccentricity of the planet might be varying with 
time. Because the time scale for orbital circularization is shorter tha n the age of the system 
by a factor of 100 (jCody fc Sasselovll2002l ; IBodenheimer et al.l 120031 ). one way to produce a 
non-zero e ccentricity for HD 209458b wou ld be to have a second planet pumping up its ec- 
centricity (IBodenheimer et al.ll200ll . 120031 ) . If th e eccentricity of HP 209 4 58b is in fact being 
pump ed u p by interactions w ith a second planet, iMiralda-Escudel (120021 ). iHolman fc Murray 
( 120051 ) and lAgol et al.l (120051 ) show that changes in the eccentricity and other orbital param- 
eters of the planet will p roduce short term oscillations in the timing of successive transits. 
Holman fc Murray! (120051 ) show that for HD 209458b the time scale for these oscillations is 
on the order of 10-100 orbital periods, depending on the orbital parameters of the second 
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planet. The radial velocity points from lLaughlin et al.l (l2005bl ) are widely spaced relative 
to the length of a single planetary orbit, and although they allow for a strong constraint on 
the average eccentricity over the five-year period covered by the data, the constraints that 
they can place on the eccentricity of an individual orbit are much weaker, as they would be 
determined by the errors on a relatively small (in most c ases less than fi ve) num ber of radial 
velocity measurements. Similarly, the observations by iDeming et al.l (j2005bl ) only span a 
single secondary eclipse and could have occurred during a period of low eccentricity. In §5.21 
we measure the locations of eight transits to check for variations in transit time that would 
indicate a time-varying eccentricity. 



Winn fc Holmanl (120051 ) propose a third possible energy source for HD 209458b. They 



argue that the planet may be in a Cassini state, in which its spin precession resonates with 
its orbital precession. In the particular Cassini state that the authors describe, the planet's 
spin axis would be tilted by close to 90° with respect to the normal of the plane of its orbit, 
and these two vectors would remain coplanar with the spin axis of the star as the system 
precessed. Because the time scale for precession is significantly longer than the orbital period 
of the planet, the direction of the planet's spin axis would be effectively fixed over the course 
of a single orbit. This would cause the substellar point to m ove from one pole to the other 
pole during a single orbital period. iBodenheimer et al.l (120031 ) calculate that a power of 10 27 
ergs s _1 is required to explain HD 209458b's enlarged radius if the planet has a dense core 
and a power of 10 26 ergs s _1 if it does not. According to lWinn fc Holmanl (120051 ). if the ratio 
of the tidal quality factor Q to the displacement Love number h is between 10 6 — 10 7 , a 90° 
obliquity would provide enough energy to explain the planet's enlarged radius. This would 
be consistent with the higher end of current estimates for the tidal quality factor of Jupiter, 
which as we noted in a previous paragraph rang e from 10 5 to 10 6 , and for the displacement 
Love number, which is estimated to be > 0.6 (jGavrilov &: Zharkovl 119771 ). Unfortunately, 
direct determination of the Cassini state of HD 209458b is beyond the capabilities of current 
techniques. 



1.2. Measuring the Radius of HD 209458b 



Since the fundamental quantity that each of these proposed models seek to explain 
is the plan etary radius, a refined v alue is of great i nteres t. The current best estimates 



come from IWittenmyer et al.l (120051 ) and IWinn et al.l (120051 ) . IWittenmyer et al.l (120051 ) use 
a combination of radial velocity data (51 published and unpublished values) and transit 
curves (Hubble Space Telescope STIS and FGS data, as well as data from several ground- 
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based observatorie s), which they fit simultaneously to obtain a radius of 1.35 ± 0.07 -RjuJE 
Winn et al.l (120051 ) also fit r adial velocity and trans it data simultaneously, but use 85 radial 
velocity measurements from Laughlin et al. (|2005a ) and transit data from the Hubble Space 
Telescope STIS observations by iBrown et al.l (120011 ) alone. The au t hors a lso include the 
constraints on the timing of the secondary transit from lDeming et al.l (j2005bl ). Fitting these 
three data se ts simultaneously, they fi nd t hat the radius of H P 209458b is 1.35 ± 0.06 -Rj up . 
The results of lWittenmyer et al.l (120051 ) and I Winn et al.l (120051 ) are virtually identical, despite 
differences in the data used in the fit. This is because the radius of the planet is determined 
by fitting the photometric transit curves (radial velocity data decreases the uncertainty in 
other parameters, and thus has a marginal effect on the b est-fit value for the radius of the 
planet, but this effect is small). The high quality of the IBrown et al.l (120011 ) HST data, 
which both authors use, relative to data from ground-based telescopes means that the HST 
data dominate any simultaneous fit of available transit curves, explaining the similarity of 
these results. 

In this work we improve on previous measurements of the planet's radius by fitting a 
new photometric data set, gathering 1066 spectra over four distinct transits with the STIS 
spectrometer on the Hubble Space Telescope. The increased wavelength range of our data 
(we use two gratings with a combined range in wavele ngth space of 290 — 1030 nm, instead of 
the single grating spanning 581.3 — 638.2 nm used by lBrown et al.ll200ll ) means that we can 
divide the spectra into ten spectrophotometric bandpasses and bin to create ten individual 
photometric t imese ries, while still maintaining a photometric precision comparable to the 
Brown et al.l (120011 ) data. The relative RMS variation for the out-of-transit data in our ten 
bandpasses ranges from 1.5 — 5.8 x 10~ 4 for a typical cadence of 40 s, with the majority of 
bandpasses closer to the lower value. This noi se is only slightly h igher than 1.1 x 10~ 4 for a 
cadence of 80 s, the typical value reported by IBrown et al.l (120011 ). 



Unlike IBrown et al.l (120011 ). who fit data spanning a single bandpass, we fit a transit 
curve over ten bandpasses simultane ously. This allow s us to break a fundamental degeneracy 
in the shape of the transit curve. As iJha et al.l (120001 ) observed, the problem arises from the 
fact that a transit is primarily described by its duration and depth. For observations in a 
single bandpass, it is possible to fit the same transit curve with a larger planet, simply by 
increasing the stellar radius and decreasing the inclination (the converse is true for smaller 
planetary radii). By observing the same transit in multiple bandpasses, we are able to deter- 
mine the inclination uniquely, independent of assumptions about the stellar and planetary 
radii. 



1 All measurements in this paper use Jupiter's equatorial radius at 1 bar, 71492 km 
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Jha et al.l (120001 ) use a similar approach in their analysis of multi-color BVRIZ photo- 
metric data fro m ground-based tele scopes, obtaining a value of 1.55±0.10 i?j up for the radius 
of HD 209458b. iDeeg et al.l (120011 ) also fit data from ground-based telescopes spanning the 
four Str0mgren bandpasses, and find a value of 1.435 ± 0.05 i?j up for the planet's radius. 
Although both authors obtain a value fo r the radius of th e planet from their fit, they differ 
in their initial assumptions. The goal of iJha et al.l (120001 ) was to obtain an improved value 
for the radius of the planet, and so they use stellar models to predict the shape of the limb- 
darkened light c urves rathe r than fitting for the limb-darkening coefficients. In contrast, the 
primary goal of IDeeg et al.l (120011 ) was to use the transit curve to probe the limb darkening 
of the primary star, and so they used a linear limb-darkening law and left the limb-darkening 
coefficient s as fr ee variables in their fit. In this work we use a similar approach to that of 
Jha et al.l (120001 ) and calculate four-parameter non-linear limb-darkening coefficients from 
theoretical models. We examine the accuracy of our limb- darkening model in more depth in 



Although we can measure the inclination independently, there is still an additional 
degeneracy between the radius of the planet (Rp), the radius of the star (-R*), and the mass 

1 /3 

of the star (M*), where Rp oc R+ oc MJ . The depth of the transit curve is a function 
of Rp/R±, so the only piece of information needed to break this degeneracy is the mass of 
the star. This can be determined by comparing the observational parameters of HD 209458, 
including its spect rum and absolute visual magnitude, with theoretical stellar-evolutionary 
models. However, ICody fc Sasselovl (120021 ) point out that the line of degeneracy between 
the stellar mass and radius from models at a constant luminosity is nearly orthogonal to 
the constraint on the stellar mass and radius from a fit of the light curves. By combining 
information from the models with the constraint from the transit curves themselves, we can 
reduce the uncertainty in the final measurement of the mass and radius of the star, which 
reduces the formal uncertainties in our measurement of the planet's radius. 

This is not an entirely new approach; early ground-based studies of HD 209458 assumed 
values for the mass and radius of the star c alculated from models, and fit for the radius o f the 
planet and inclination of its orbit alone (jCharbonneau et al.l |2000| ; iHenrv et al.l 120001) . A s 
higher-quali t y dat a fro m the Hubble Spac e Teles cope became available, iBrown et al.l (120011 ). 
Winn et al.l (120051 ) and IWittenmver et al. (hoosh assume a value for the mass of the star and 



use this data (originally published by lBrown et al.ll200ll ) to measure the value for the stellar 
radius directly from the transit. The quality of the data is important because it is the shape 
of the ingress and egress that breaks the degeneracy between inclination and stellar radius 
for data in a single bandpass. Although our data allow us to fit for the stellar r adius directly 



(given a value for the stellar mass) , with an accuracy greater than that of IBrown et al. 



( 120011 ) . in our analysis we incorporate the additional joint constraint on the stellar mass 
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and radius from models (jCody fc Sasselovll2002l ) in order to minimize the uncertainty in the 
measurement of the planetary radius. 



2. Observations and Data Reduction 

We obtained twenty HST orbits of STIS spectra, grouped into four "visits" of five 
consecutive orbits each, during which the telescope was pointed continuously at the star. 
Each visit was timed so as to span one complete transit of HD 209458b. These observations 
were obtained during the period from UT 2003 May 3 to July 6 (see Figure [Q as part of 
the GO-9447 program. For a more detailed discuss ion of the sui t abilit y and performance of 



the STIS spectrometer for transit photometry see iBrown et al.l ( 120011 ). We observed using 
both the G430L and G750L gratings, with two visits in each grating. Together these two 
filters cover a combined range of 290 — 1030 nm, with some overlap between filters around 
530 nm (see Figure [2]) • In our analysis we sub-divide the spectrum and bin it into five equally 
spaced bandpasses in wavelength space within each grating, yielding a total of ten distinct 
bandpasses (see Figure [3]). We note that there is significant fringing at wavelengths longer 
than 800 nm from the internal interference of the thin STIS CCD, but we find this fringing 
has a negligible effect on the binned data. We gathered spectra with a 22 s integration time 
at a cadence of 42 s in the G430L grating and a 19 s integration time and 39 s cadence in the 
G750L grating within the part of each spacecraft orbit where the star was visible, yielding 
excellent time resolution during the crucial ingress and egress periods. 

We assigned wavelengths to our data by cross- correlating our spectra (using the wave- 
length solution returned by the automated STIS pipeline) with the model spectrum of 
HD 209458 described in SJ3J We performed this analysis over 20 individual bandpasses 
spanning the spectral range of the data. We found that the typical offset did not vary by 
more than one pixel in time or with spectral region, and hence we calculated the mean offset 
value and shifted the STIS wavelength solution accordingly. 

We create our timeseries for each bandpass by binning the portion of each individual 
1024 x 64 pixel spectrum within the desired range of wavelengths to create a single photo- 
metric measurement. We optimized the size of our bin in the cross-dispersion direction to 
minimize the contribution from read noise and scattered light distributed across the detector 
array. This was particularly important in the regions of the spectrum where the signal was 
small, as the additional signal obtained from a larger bin in the cross-dispersion direction 
was minimal compared to the increase in noise. 



We optimized our aperture sizes in the cross-dispersion direction for 50-pixel-wide sec- 
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tions spanning the range of the spectrum. In order to determine the optimal size, we picked 
an aperture, binned the data to create a timeseries, and measured the relative variation in 
the out-of-transit data for that timeseries. We repeated this process systematically for all 
possible aperture sizes, ranging from 3 — 64 pixels, and selected the aperture that minimized 
this variation. The optimal aperture sizes in the cross-dispersion direction for each of our 
segments ranged from 5 — 31 pixels, with a median size of 19 pixels. The narrowest apertures 
were in the low-flux regions at the edges of the spectra and the wider apertures in the central 
high-flux regions. 

Once we had determined our optimal apertures and binned our spectra within each 
bandpass to create a photometric timeseries, we found that there were several prominent 
trends remaining in the data. First, the flux values increased gradually over time during 
each visit, so that the average flux before the beginning of a typical transit was 0.1% lower 
than the average flux after the transit. This trend was particularly pronounced in data taken 
during the first spacecraft orbit of each visit, which was sometimes as much as 0.3% lower 
than data at the end of the visit. We attributed this trend to thermal settling of the telescope 
at its new pointing. Second, the first point of each spacecraft orbit was typically 0.2% lower 
than subsequent points. Third, the data also showed a trend within individual spacecraft 
orbits, with a typical range of 0.1 — 0.2% between minimum and maximum values, that was 
repeated consistently from one orbit to the next. Although it is not clear what caused these 
variations, they were found to correlate with spacecraft orbital phase. This indicates that 
they were caused by changes in the instrument and not intrinsic variations in the star. 

In order to remove these trends, we used two steps. First, we chose to discard the first of 
the five orbits in each visit and the first point in each of the remaining orbits. The first orbit 
and initial points of subsequent orbits consistently exhibited the largest systematic effects, 
and discarding these points still left sufficient data before and after the trans it to measure 



an accurate value for the out-of-transit flux from the star. iBrown et al.l (120011 ) removed the 
same points from their data, which showed similar trends. In the second step, we removed 
the remaining trends by fitting the out of transit points in each transit simultaneously with 
a third-order polynomial function of spacecraft orbital phase and a linear function of time, 
and dividing the data by the resultant function. This fit was done individually for each of 
the four visits, in each bandpass for that visit. 

The normalized timeseries for each transit (Table [1]) is quite flat, although as we will 
discuss in §5.2[ there are still some remaining systematic effects visible in the combined plots 
(the plot for each bandpass includes data from two transits overplotted) in Figure [3J The 
RMS variation in relative flux for the out-of-transit data in each of these bandpasses ranges 
from 1.5 — 5.8 x 10~ 4 , with the noisiest bandpasses at the low-flux ends of the spectrum. The 
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time coverage of the data for the combined transit curve is quite good, as shown in Figure 

m 

When we compared the RMS scatter of the out of transit data to the expected photon 
noise for the data in each bandpass, we found that the RMS scatter was on average 21% 
higher than the photon noise. We examined several potential noise sources. The read noise 
for the STIS detector is 7.75 e~ pixT 1 RMS at a gain of four, with 3800 pixels in an average 
bandpass. A typical bandpass has a flux of 10 8 e~ and a corresponding Poisson noise of 
10 4 e~ RMS, so read noise would only decrease the SNR by 0.1%. There is also a constant 
background contribution from a combination of light reflected from the Earth, zodiacal light, 
and geocoronal emission. This background contributes up to 0.06 e~ s~ 4 pix -1 , for a total 
contribution of 10 4 e~. This is only 10~ 4 of the total signal, so its contribution to the Poisson 
noise is negligible. Although the magnitude of this background may vary over time as the 
position of the telescope changes relative to the various sources, this is unlikely to be the 
source of the additional variation, for two reasons. First, although some components of the 
systematic variations in the data correlate well with changes in the width and location of the 
point spread function of the spectrum in the cross-dispersion direction, they do not correlate 
with the relative positions of the Sun and moon, two likely sources of scattered light within 
the instrument. Second, when the data is binned using increasingly narrow slits in the cross- 
dispersion direction, it is clear that most of the systematic variations are contributed by the 
central 20 pixels around the peak of the point spread function, making scattered light in the 
instrument an unlikely explanation. The most plausible origin is the change in the focus 
of the instrument and location of the source in the slit, which do not precisely repeat from 
orbit to orbit. In order to properly account for all noise sources in our error estimates for 
the points in our timeseries, we chose to set the error for all points within a given bandpass 
as the standard deviation of the out of transit points for the data in that bandpass. 



3. Fitting the Transit Curve 



We use the complete analytic formula given in lMandel fc Ago 



(2002), without 



approx- 



imations, to calculate our transit curves. The expression given in iMandel fc Agoll (120021 ) is 
a function of six dimensionless parameters, including the ratio of the planetary radius to 
the stellar radius, the impact parameter in units of stellar radii, and four nonlinear limb- 
darkening coefficients. We would like to fit for the first two variables, as these variables 
determine the best-fit values for the mass and radius of the star, the radius of the planet, 
and the inclination of the planet's orbit relative to the observer. We have designed our fit- 
ting routine to use the latter four parameters as its input, calculating the two dimensionless 
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parameters for the analytic expression from iMandel fe Agoll (120021 ) . 



Rather than fitting for the period P and ini tial transit time Tg sim ultaneously, we fix 
these p arameters to an initial estimate, taken from lBrown et al.1 (120011 ) and lCharbonneau et al 
(120031 ). and fit for the best planet radius, inclination, stellar mass, and stellar radius. We 
then go on to find P and Tg using the method described in §H We then take the new best-fit 
P and Tq and repeat our fits for planet radius, inclination, stellar mass, and stellar radius. 
We iterate this process until our best-fit values for all six parameters converge to a consistent 
solution. Because we are using a relatively simple fitting routine (downhill simplex) to avoid 
the need for derivatives of our transit function, this iterative fitting ensures that the function 
converges to the true global minimum. We assume the orbital eccentricity is zero, for the 
reasons discussed in SJTJ 

As discussed in the introduction, we must make an initial assumption about either the 
mass or the radius of the star in order to break the degeneracy between these parameters 
and the radius of the planet when we fit our tr ansit curve. We constrain the mass and radius 
of the star based on ICody fc Sasselovl (120021 ). where R pre d is the radius predicted by the 
theoretical mass-radius relation discussed in SJT] for a given stellar mass: 



R 



pred. 



1.18 



- 1.06 



R & 0.96 \M Q 

We implement this constraint by adding two terms to the chi-squared function where we treat 
M* and R* as Gaussian random variables with standard deviations equal to the published 
errors. The inclination and planetary radius have no such constraints. The goodness-of-fit 
parameter is given by: 




(2) 



where rrii is the i th measured value for the flux from the star (with the out-of-transit points 
normalized to one), pi is the predicted value for the flux from the theoretical transit curve, 
and a m is the error on each individual flux measurement, which we take to be the standard 
deviation of the out of transit points for that bandpass, as discussed in §2J M* and R+ are 
the fitted values for the mass and radius of the star, and we constrain the mass of the star 



using the value from ICody fc Sasselovl (120021 ). 1.06 ± 0.13 M w . The er ror on the predicted 
radius of the star, ±0.055 R & , also comes from lCody fc Sasselovl (120021 ) . 



When calculating our transit curves, we use the nonlinear limb-darkening law defined 
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in Oaretl (|2000|): 



J(/i) = !-$>(!- 



n/2\ 



n=l 



where 



cost? 



(3) 



(4) 



In this case, 9 is defined as the angle between the normal from the surface of the star and the 
direction of the observer. We derive our four-parameter nonlinear limb-darkening coefficients 
from theoretical models calculated by R. KuruczB (personal correspondence, August 2005) 
using an effective temperatu re of 6100 K, log(g) = 4.38, and [M/H] = +0.014 for the 
star (iValenti fc Fischer! 120051 ). The model spectra used in these calculations include both 
continuum processes and line absorption fcuruczl [2005l l The limb-darkening coefficients 
for each bandpass are determined by taking the weighted average of the limb-darkening 
coefficients at individual wavelengths, with the weights determined by the measured flux at 
that wavelength (see Tableland Figure[2]). This includes the wavelength-dependent response 
function of the STIS instrument. As expected, we find that our limb-darkened curves are 
similar to those of the Sun at longer wavelengths, but exhibit greater center-to-limb variation 
at short wavelengths where the relative difference in intensity between blackbody curves at 
the temperature of the two stars is greatest (see Tableland Figure [3]). 

To determine the best-fit radius for the planet, we evaluate the \ 2 function over all 
ten bandpasses simultaneously, using the same values for the planetary radius, stellar mass 
and radius, and inclination. Minimizing this function, we find that the best-fit value for the 
radius of the planet is 1.320 ± 0.025 i?j up , with an orbital inclination of 86?929 ± 0?010. The 
best-fit mass and radius for the star are 1.101±0.066 M and 1.125±0.023 R Q , respectively. 
The reduced \ 2 f° r this minimized function is 1.28. 

As a check, we also fit for the limb-darkening coefficients directly using the data. Because 
we would like to compare the shapes of these curves to the curves derived using our model 
limb-darkening coefficients, we set the stellar mass equal to its best-fit value from Table H] 
and continued to include an additional term in the x 2 function for the mass-radius relation 
from lCody fc Sasselovl (120021 ) . We also chose to fit for the limb-darkening coefficients using 
a quadratic limb-darkening law iMandel &: Agoll (120021 ) instead of the four-parameter non- 
linear law we used earlier, as the four-parameter coefficients are degenerate at the precision of 
our data. In order to fit for all 24 variables, including 20 limb darkening coefficients, reliably 
we use a Markoy Chai n Monte Carlo code d eveloped following the methods described by 
Ford et al.l (120051 ) and iTegmark et al.l (12004 ). and check our results using our original x 2 



2 Available at |http://kurucz.harvar cl.edu/stars/hd2 09458 
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minimization routine. When we include the uncertainty in the stellar mass from Table HI 
we find R+ = 1.137 ± 0.024, R P = 1.342 ± 0.032, and i = 86?75 ± 0?14. As expected, this 
new fit is able to remove some of the trends in the residuals, with the largest remaining 
discrepancies around the periods of ingress and egress (see Figure H]). This is not surprising, 
as the uncertainties in the limb-darkening laws predicted by models are greatest at the edges 
of the star. 

Although these new values are a better fit to the data, the uncertainties in the best-fit 
parameters are also correspondingly larger. This new fit is also more sensitive to residual 
noise correlated with spacecraft orbital phase (discussed at the end of §2j) , which can alter the 
shape of the ingress and egress. As a result, we prefer the more robust fit with model-derived 



l imb-d arkening; coefficients for subsequent analysis. In an independent study, iTingley et al. 



(120061 ) also compared the observed limb-darkening in published transits for HD 209458 to 
the predicted limb-darkening from the ATLAS models, and concluded that although there 
were some discrepancies between the two curves these discrepancies had a negligible effect 
on the measurement of the system parameters. We find that in this case our values for the 
radius of the planet differ by less than la between the two methods. 



4. Fitting the Ephemeris 



In order to measure an accurate value for the planetary ephemeris, we included four 
normalized transit curves^ derived from STIS data taken between UT 2000 April 25 and May 
13 in our analysis. This data was taken using a higher-resolution grating (G750M instead 
of G430L and G750L) and spans the wav elength range from 581.3 — 638.2 nm, instead of 
the 290 — 1030 nm spanned by our data. iBrown et al.l (120011 ) binned this data in a single 
ban dpass over the enti re available wavelength range. As discussed in §1.21 the precision of 
the IBrown et al.l (120011 ) transit curves and our transit curves are comparable. The authors 
observe the same kinds of systematic effects in their data as we do in ours, although the 
orbit-to-orbit residuals in our data are larger, and they normalize their data with a linear 
function of time from first observation and a fourth-order polynomial function of orbital 
phase to remove th ese effects. We use these normalized transit curves, taken directly from 
Brown et al.l (120011 ). in our fits. 



Rather than using the quadratic limb-darkening coefficients published by IBrown et al. 



(2001), we derive new four-parameter nonlinear limb- darkening coefficients for the lBrown et al. 
(120011 ) data using the same method as described in Sj3l As before, we define our limb- 



3 Available at 



http://cfa-www.harvard.edu/~dcharbon/frames.html 
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darkening coefficients for the bandpass as the weighted average of the nonlinear limb- 
darkening coefficients at each wavelength, with the flux from a typical spectrum at that 
wavelength as the weight. Although all four transits were observed with the same grating, 
there was a database error in the location of the subarray for the first transit, and as a result 
the red part of the spectrum was not entirely contained on the subarray. We account for this 
misalignment by calculating one set of limb-darkening coefficients for the first visit and a 
different set for the other three visits. For the first visit our four-parameter non-linear limb- 
darkening coefficients are c\ = 0.5838, c 2 = —0.3430, c 3 = 0.9183, and c 4 = —0.4240. For the 
subsequent three visits c x = 0.5845, c 2 = -0.3443, c 3 = 0.9179, and c 4 = -0.4236. As be- 
fore, we check the standard linear wavelength solution from the HST pipeline by comparing 
these spectra to model spectra for HD 209458. Using the sodium lines at 588.995 nm and 
589.592 nm as our calibration we find the HST wavelength solutions for this line are shifted 
by 1 pixel (0.06 nm) relative to the rest wavelength of the line, and correct the wavelength 
solutions accordingly before calculating our limb-darkening coefficients. 

We fit for the locations of each of the eight transit centers individually using a standard 
X 2 function and the best-fit values for planetary radius, inclination, and stellar mass and 
radius from §3J We bin the data from 2000 into a single bandpass for each transit, and 
leave the data from 2003 in the same bandpasses used in the previous analysis, fitting 
simultaneously over all bandpasses. After deriving the best-fit locations for each transit 
center, we fit for the period and initial transit location Tc by plotting the number of transits 
versus the location of each transit in HJD and fitting the points with a line. The constant 
coefficient of this linear fit gives us the best-fit Tc, while the slope of the line gives the period. 
We then repeat our earlier fits for Rp, i, i?*, and M* using this new value for the period and 
initial transit time, iterating until all values converge to a consistent result. 



Discussion 



5.1. System Parameters 



The best-fit values for the planet's radius and orbital inclination are listed in Table 
IH Figure [3] shows the data in each bandpass with best-fit transit curves overplotted, and 
Figure H] shows the residuals from the fit using the theoretical four-parameter non-linear limb- 
darkening coefficients. The new value for Rp, 1.320 ± 0.0 25 -R,t„ p , has uncert ainties that are 
less than half that quoted in IWittenmyer et al.l (120051 ) and I Winn et al.l (120051 ) . Although our 
radius is slightly smaller than the values given by these two authors, it is entirely consistent 
within their uncertainties. Using the mass given in lLaughlin et al.l (l2005bl ). 0.69 ±0.06 Mj up , 
and scaling it appropriately to account for the different value for the stellar mass from our 
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fits, we calculate Mp = 0.64±0.06 Mj up . This gives a density for the planet of 0.26±0.04 pj up 
or 0.345 ±0.05 g/cm 3 . 

We find that the inclination of H P 209458b's orbit is 86?929 ± 0?010. Although our 



for the inclination given by IWittenmyer et al 



inclination is outside the er rors given by I Winn et a l. (200 51) (there are n o independent errors 



2005 



Winn et al.l (120051 ) note that the incli- 



nation depends on their particular choice of limb-darkening law. In this case, we find that 
we measure a very small error for the inclination when we assume the limb-darkening coef- 
ficients are known from models. Unfortunately, much of this improvement is lost when we 
allow the limb- darkening coefficients to vary in the fit. In this case, our use of ten individual 
bandpasses is no longer an advantage, as limb-darkening adds two degrees of freedom to our 
fit for every bandpass we include. This makes it possible to fit the data with a significantly 
larger range of inclinations. It also makes the fit more sensitive to correlated noise in the 
data, as discussed at the end of $2j For this reason, we prefer the more robust fit using 
model limb-darkening coefficients, which provide a reasonably good fit to the data and do 
not have a significant effect on the value we measure for the radius of the planet. 

We note that the value for the radius of the planet that we obtain from our fits di- 
rectly depends on our initial assumptions about the mass and radiu s of the sta r. Al- 



t ho ugh we use the same init ial estimate for the mass of the planet as IWinn et al.l (120051 ) 
and IWittenmyer et al. J200A we include an additional restriction on the stellar radius as a 
function of the stellar mass (jCody &: Sasselovll2002l ). We find that best-fit values for the mass 
and radius of the star from our fit are 1.101 ± 0.066 and 1.125 ± 0.023 -R^, respectively. 



These values differ from the best-fit values cited in ICody fc Sasselovl (120021 ) that we used in 
our initial constraints because the transit curve provides additional information about the 
stellar radius as a function of stellar mass, as discussed in 81.21 T hese values are within 
the quoted errors for the measurements from I Cody fc Sasselovl (120021 ). and are also virtually 
ident ical to the most recent in dependent measurements of the mass and radius of HD 209458 
from IValenti fc Fischer! (120051 ). who find a stellar mass and radius of 1.11 ±0.18 M and 
1.122 ±0.055 R Q . 

We also repeat our fits without any constraints on the stellar radius. In this case, the 
best-fit stellar radius is de termined by the value we assume for the mass of the star. Using the 
relatively small mass from I Cody fc Sasselo v (2002 ), we find that .R * = 1.11 ±0.04 R Q , a value 
that is smaller than the one measured by ICody fc Sasselovl (120021 ) but still within 2a of that 
measurement. The radius of the planet in this fit is 1.31 ± 0.05 -Rj up . The increased errors 
indicate that it is the joint constraint on the mass and radius of the star which produced 
the most signification reducti on in the uncertai n ty in the planet's radius in our earlier fits. 
When we use the mass from IValenti fc Fischerl (120051 ). we find R* = 1.13 ± 0.06 R Q , and 
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R P = 1.33 ±0.07 R Jup . 

In our analysis so far we have assumed that there was a single universal value for 
Rp, R+, M*, and i. However, R P could vary from bandpass to bandpass, as the presence 
of absorption lines in the planet's atmosphere might increase the depth of the eclipse in 
particular bandpasses. This is particularly relevant for the long-wavelength bands in the red 
grating, where model at mospheres for HD 209458b predict some of the strongest absorption 



bands from water vapor (jBrownll200ll ). We check the validity of our previous assumption by 



fitting each bandpass individually. 

Our initial results were surprising; we measured a variation in the best-fit values for 
Rp on the order of ±0.02 i?j up . Although it is tempting to dismiss this variation as within 
the la errors for Rp, it is important to remember that these errors are dominated by the 
error in the values of M+ and R+, quantities that do not vary from bandpass to bandpass. 
Because we applied the same constraints in each of our fits over individual bandpasses, we 
would expect the relative value of Rp within each bandpass to remain consistent over much 
smaller scales than the overall error for the measurement. 

Closer examination of the best-fit values for the other variables in the fit revealed that 
each value of Rp was associated with a different value for the inclination. Of course, incli- 
nation cannot vary with bandpass, and so we repeated the fits while fixing the inclination to 
the best-fit value from the simultaneous fit of all ten bandpasses. Because we are interested 
in the relative uncertainty in the depth of the transit between bandpasses, we also fix the 
mass and radius of the star to their best-fit values. As a result, the best-fit radii in Table 
[5] reflect the relative uncertainty in the measurement of the radius in a particular bandpass, 
and not the total uncertainty in the radius measurement. This decreased the variations in 
i?P between bandpasses significantly, to an average of ±0.003 -Rj up or 210 km, with the 
largest shifts in the 293 — 347 nm and 922 — 1019 nm bandpasses. 

As we noted earlier, some of the residuals in the data come from the discrepancy be- 
tween our model limb-darkening coefficients and the observed limb-darkening of the star. In 
order to check the effect that the limb-darkening coefficients might have on our results, we 
repeated our fits using the best-fit quadratic limb-darkening coefficients from the Markov 
Chain Monte Carlo fit described in |J3J We found that although the change in radius for a 
particular bandpass was not always the same as it was for the fits using model limb darkening 
coefficients, the overall level of variation between bandpasses was virtually identical. 

Although we expect absorption lines to produce some differences between bandpasses, 
a change of 0.003 i?j up in the effective radius of the planet would correspond to an average 
change of 0.07% in the relative depth of the transit over the entire 50 — 100 nm wide bandpass. 
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By comparison, the sod i um ab sorption line in HD 209458b's atmosphere at 589 nm discovered 
by ICharbonneau et al.1 (120021 ) only produces a change of 0.02% in the relative depth of the 
transit over a range of 1 nm in wavelength space. Although molecular absorption bands could 
produce a significant change in the depth of a transit over a wider range in wavelength space, 
the lines would have to be more than three times stronger than the sodium line discovered by 
Charbonneau et al.l (120021 ) and cover most of the bandp ass to explain the observed variations . 



Based on model tran smission spectra for HD 209458b (ISeager fc Sasselovll2000l ; lBrownll2001 



Hubbard et al.ll200ll ). there are only a few metal absorption lines predicted in most of our 
bandpasses, and the only absorption bands are from water absorption and are located in the 
longest-wavelength bandpasses of our data. If these models are accurate, absorption at the 
level of the variations we measure is unlikely. 

We can also perform a more general estimate of the potential strengt h of absorption 



l ines. U sing our new value for the radius of the planet and the mass given by lLaughlin et al. 



(l2005bl ) (appropriately scaled to reflect our smaller va l ue for the stellar mass), and the 
brightness temperature at 24 fim given by iDeming et al.l (j2005bl ) , we calculate a scale height 
of 450 km, or 0.006 -Rj up . If the opacity in absorption lines was effectively one over a range of 
one scale height above the average radius we measure, we would expect to see an increase of 
0.006 Rj up in the best-fit radius in the region of the absorption line. This is twice the average 
variation we observe. This means that it is possible, if unlikely g i ven w hat we already know 
from the detection of s odium absorption by ICharbonneau et al.l (120021 ) and detection limits 
given by iNarita et al.l (120051 ) in this wavelength range, that this level of variation could be 
the result of absorption in the atmosphere of HD 209458b. However, it is entirely possible 



that s maller variations, on the level of the Na absorption detected by ICharbonneau et al. 



(120021 ). might be detectible in the data. We plan to present the results of a thorough search 
for spectroscopic transmission features in an upcoming paper. 



5.2. Orbital Ephemeris 



We derive a period of 3.52474859 ±0.00000038 days (±0.033 s) for the planet and initial 
transit epoch T c = 2452826.6285 2 1 ± .000087 HJD (±7.5 s). This period is consistent 
with the period fromlBrown et al.1 (1200 ll ). 3.52480 ± 0.00004. We compare our Tc with the 



results of iBrown et al.l (l200ll ). repeating our fit using the first transit in our dataset as our 
initial transi t time and find T c = 2451659.936739 ± 0.000087 HJD. This differs from the 
Brown et al.1 (120011 ) value by 1.1 x 10 -5 days or 1 s, and was also well within the published 



error of 10 4 days (see Table [7]). 



Our period and time of center of transit are entirely consistent with the values from 
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Brown et all (I20011). which use HST STI S data alone. However, our period differs from the 
period g iven by IWittenmver et all (120051 ) by 26 s, or 17a. When we use the period and T c 



Wittenmver et all too£ to calculate the predicted times for the UT 2003 July 5 



given m 

transit, the value differs by 118 s or 10cx from our best-fit time. We note that our data is more 
recent and has better time coverage within each transit than the HST Fin e Guidance Sensor 



data from 2001 and 2002 that dominates the period measurement by IWittenmver et al. 



(2005). 



Although the period for the planet is consistent with previous results by iBrown et al. 



( 120011 ). we initially found significant discrepancies between the predicted and measured tran- 
sit times for individual transits (see Table [6]) . We derive the errors for our timing measure- 
ments using a bootstrap Monte Carlo method, in which we randomly drew measurements 
with replacement from the data set of a given transit to simulate a data set with the same 
number of points and noise properties as the original, we then fit for the time of center of 
transit. We repeat this process until we obtain an approximately Gaussian distribution of 
best-fit transit times, and take the standard deviation of this distribution as the error on 
our timing measurement. Using this method, we find errors ranging from 3 — 6 s. When we 
calculate a new P and T c using these errors, we find that our best-fit transit times differ 
from their predicted values by as much as 6a. However, as Figure [5] demonstrates, there are 
still some residual trends in the data that are not taken into account by this method. 

We examined the flux residuals from the transit light curves and concluded that there 
were residual trends in the data that might alter the measured transit times. Although some 
of these trends are the result of our choice of limb-darkening coefficients, as discussed in §3j 
there are also residuals from our procedure to correct for the instrument variations (described 
in Sj2]) that might affect the results. Because systematic shifts over an entire spacecraft orbit 
will have the largest effect on the transit times, we measured the typical offset of each orbit 
in the residuals plotted in Figure HI We found that the typical offset for an individual orbit 
was 0.0001 in relative flux. In order to quantify the effect that this kind of orbit-to-orbit 
variations in flux might have on the timing measurements, we created artificial data sets 
for each of the eight transits and shifted the points in each individual spacecraft orbit by 
a Gaussian random variable with a standard deviation of 0.0001 in relative flux. For this 
level of variation, we calculated a standard deviation of ±14s in the best-fit transit times. 
Adopting this value as the error on our timing measurements, we find that all shifts are 
within 2a of their predicted values (see Figure EJ). The errors for P and Tc are derived 
assuming a constant error of ±14s on each individual transit time. 
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Conclusions 



We used multiple-bandpass photometry, the stellar mass-radius relation from lCody fc Sasselov 



(120021 ). and theoretical models for limb darkening to significantly improve the estimates of 
the radius and orbital inclination of HD 209458b. We find that the radius of HD 209458b 
is 1.320 ± 0.025 -Rj up , a fact or of two more precise than previous measurements. Using the 
mass of the planet given by lLaughlin et al.l (j2005bl ) and scaling it appropriately to account 
for our smaller value for the stellar mass, we find a density for the planet of 0.26 ± 0.04 pj up 
or 0.345 ± 0.05 g/cm 3 . The planet's inclination is 86?929 ± 0?010, a factor of three more 
precise than previous measurements. 

We also improve upon previous measurements of the planet's orbital period P and 
time of center of transit Tq- Using a combination of STIS data spanning four transits in 
2000 and four transits in 2003, we find that P = 3.52474859 ± 0.00000038 days (±0.033 s) 
and T c = 2452826.628514 ± 0.000087 HJD (±7.5 s). Although the quoted uncertainties in 
the period for HD 209458b derived by IWittenmyer et all (120051 ) are smaller, our data has 
better signal-to-noise and time coverage than the Hubble Space Telescope FGS data which 
dominates their me asurement. Additi o nally, the large (10er) discrepancy between the transit 
times predicted by IWittenmyer et al.l (120051 ) and the transit times we measure for the 2003 
data ind icate that there may be systematics that are not taken into account in the error 
given by IWittenmyer et all (120051 ) . 



When we compare the observed timing of individual transits to the predicted values, 
we see no discrepancies larger than 2a. Our flux residuals a nd the errors on ou r timing 
measurements are modestly larger than the errors quoted by iBrown et all (120011 ) . but we 
extend the time frame of these results by three years and double the number of transits 
observ ed at a comparable level of precision. Using the equations from iHolman fc Murray 
( 120051 ) . and setting a limit of 42 s or 3<j on potential timing variations, we are able to 
place a limit of 4 Mj up on the mass of a potential second planet in this system, where we 
have assumed either a relatively short period and circular orbit, or a longer period (up to 
100 days) and larger eccentricity ( up to e = 0.7). We note that this limit is consistent with 
the smaller limit of 0.3 Mj up that lLaughlin et all (l2005bl ) derive from radial velocity data 
for this system. 

Given the success of this method, and the need for accurate characterization of extra- 
solar planet radii, it would be useful to apply this method to other transiting hot Jupiters. 
Despite the failure of STIS, recent results by T. Brown (personal communication, 2005) 
have demonstrated that spectra with a cadence and signal-to-noise comparable to that of 
STIS can be obtained using the grism spectrometer on HST's Advanced Camera for Sur- 
veys. This opens up the possibility of similar multiple-bandpass studies for other transiting 
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planets, an avenue that we plan to pursue in the future. There are currently five known 
extrasolar planets around stars bright e nough for this kind of study: HP 209458b, TrES-1 



200 



6; 



2006 ri). HP 189733b flBouchv et al 



(lAlonso et all 12004 ; ISozzetti et al.ll2004r). HP 149026b flSato et al.ll2005l : ICharbonneau et al. 



20051 iBakoset al.l bood ). XO-lb JMcCullough et al. 



Holman et al. 20061 ). and it is likely that several more will be uncovered in the coming 



years by a combination of wide-field transit surveys and quick-look radial velocity surveys. 

We are grateful to R. Kurucz for his assistance in calculating the limb-darkening for 
HP 209458. Support for program number HST-GO-09447 was provided by NASA through a 
grant from the Space Telescope Science Institute, which is operated by AURA under NASA 
contract NAS5-26555. H.K. was supported by a National Science Foundation Graduate 
Research Fellowship. 
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Fig. 1. — Normalized timeseries spanning the entire wavelength range of the grating for 
the four transits observed. Transits occurred (from top to bottom) on UT 2003 May 3 
(HJD 2452763.7, G430L grating), May 31 (HJD 2452791.4, G750L grating), June 25 (HJD 
2452816.1, G430L grating), and July 5 (HJD 2452826.6, G750L grating). Each successive 
transit is offset by 0.008 in relative flux. 

Table 1. Photometry of HD 209458 



^center (mil) HJD Relative Flux Uncertainty 



320.1 


2452763.055640 


1.000531 


0.000609 


320.1 


2452763.056126 


1.000054 


0.000609 


320.1 


2452763.056612 


1.000691 


0.000609 


320.1 


2452763.057098 


1.000254 


0.000609 


320.1 


2452763.057585 


1.000361 


0.000609 



Note. — The quoted uncertainties are the standard deviation of 
the out of transit data in each bandpass. We intend for this Table 
to appear in entirety in the electronic version of the Astrophysical 
Journal. A portion is shown here to illustrate its format. The data 
are also available in digital form from the authors upon request. 
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Fig. 2. — TOP: Two typical extracted one- dimensional spectra. The G430L spectrum spans 
the wavelength range from 293 — 567 nm, and the G750L spectrum spans the wavelength 
range from 532 — 1019 nm. The vertical lines denote individual bandpasses within each 
grating, with some overlap between gratings around 550 nm. Fringing from the detector is 
significant for wavelengths longer than 750 nm. BOTTOM: Limb-darkened light curves for 
the ten bandpasses described above, as discussed in $3j plotted using the best-fit parameters 
from Table H] and model four-parameter non-linear limb- darkening coefficients. 
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Fig. 3. — Normalized data for the ten bandpasses shown in Figure [2], with theoretical transit 
curves using the best-fit parameters from a simultaneous fit of all bandpasses (Tables H] and 
[7]) and model four-parameter non-linear limb- darkening coefficients overplotted. Note that 
each bandpass contains data from two separate visits, consisting of four spacecraft orbits 
each (visits are plotted individually in Figured] for reference). Each successive transit curve 
is offset by 0.004. Although the transit curves are a good fit for most of the data, there are 
some systematic deviations on timescales comparable to that of a HST orbit (see Figure Hj) 
that are discussed in §5.2! 
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Fig. 4. — These are the residuals in each bandpass for our best-fit parameters. Each suc- 
cessive bandpass is offset by 0.0028. As noted in §5.2[ the relative fluxes shift by 10 -4 on 
average from one orbit to the next. Each continuous section of data consists of two spacecraft 
orbits from different visits. The five redmost (upper) curves were all gathered simultane- 
ously, resulting in the correlated variations among the five bandpasses. Similarly, the five 
bluemost (lower) curves were gathered simultaneously. The solid curves overplotted are the 
difference between transit curves calculated using fitted quadratic limb- darkening coefficients 
and transit curves using the nonlinear limb- darkening coefficients derived from the model by 
Kurucz. 
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Fig. 5. — These are the O-C residuals for all eight transits using the period and Tc derived 
in this work. The errors shown are based on an extrapolation of the systematic spacecraft 
orbit-to-orbit variations in the data; see §5.2! for a complete description. The dashed lines 
are calculated from the uncertainties in the measurements of P and Tc- The two lower plots 
show the relevant regions of the upper plot in more detail. 
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Table 2. Nonlinear limb-darkening coefficients from model 



A (nm) ci C2 C3 C4 



293 - 


347 


-0.1015 


0.5547 


0.6096 


-0.2814 


348 - 


402 


0.0284 


0.4248 


0.6646 


-0.3450 


403 - 


457 


0.2022 


0.1101 


0.9690 


-0.4801 


458 - 


512 


0.3765 


-0.0119 


0.8863 


-0.4504 


512 - 


567 


0.4957 


-0.2057 


0.9157 


-0.4340 


532 - 


629 


0.5566 


-0.2972 


0.9190 


-0.4288 


629 - 


726 


0.6239 


-0.4176 


0.8889 


-0.4026 


727- 


824 


0.6495 


-0.4916 


0.8722 


-0.3844 


825 - 


922 


0.6623 


-0.5338 


0.8411 


-0.3661 


922 - 


1019 


0.6535 


-0.5323 


0.8142 


-0.3566 



Table 3. Comparison of limb-to-central intensity 



A (nm) 


Our model 


Best Fit 


Solar 


293 - 347 


0.219 


0.122 


0.066 


348 - 402 


0.227 


0.172 


0.083 


403 - 457 


0.199 


0.226 


0.115 


458 - 512 


0.199 


0.314 


0.152 


512 - 567 


0.228 


0.357 


0.188 


532 - 629 


0.250 


0.317 


0.216 


629 - 726 


0.307 


0.468 


0.272 


727 - 824 


0.354 


0.521 


0.317 


825 - 922 


0.397 


0.524 


0.363 


922 - 1019 


0.421 


0.479 


0.393 



Note. — The second column lists the limb-to- 
central intensity ratio for the four-parameter non- 
linear limb-darkening coefficients that we use in 
our fits, calculated from a model for HD 209458. 
The third column lists the intensity ratio for the 
quadratic limb-darkening coefficients that we ob- 
tain by fitting the data directly, and the fourth 
column lists the intensity ratio for theoretical four- 
parameter non-linear lim b-darkening coefficients 
for the Sun, calculated by I Claret! j2000ll . 
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Table 4. Comparison between best-fit values and results from previous works 



Study R P (Rjup) 


Inclination (°) 


M* (M ) 


R* (R©) 


Wittenmver et al. f2005~) a 1.35 ± 0.07 
Winn et al. f2005) b 1.35 ± 0.06 
This Work a 1.320±g;g|| 


86.668 
86.55 ±0.03 
86.929I° ;» 


1.09 ± 0.09 
1.06 ± 0.13 
1 101+ 066 


1.15 ±0.06 

1 icr+0.05 
1 - lo -0.06 
i 19 c-+0.020 
l.izo_ Q23 



a Used stellar mass-radius relation from lCodv fc Sasselovl J2002h 
b Assumed value for the stellar mass from lCodv fc Sasselovl <2002h 



Table 5. Best-fit Planet Radius in Individual Bandpasses 



A (nm) 


Radius (-Rj up ) 


293 - 


347 


1.3263 ±0.0018 


348 - 


402 


1.3254 ±0.0010 


403 - 


457 


1.3200 ±0.0006 


458 - 


512 


1.3179 ±0.0006 


512 - 


567 


1.3177 ±0.0010 


532 - 


629 


1.3246 ±0.0006 


629 - 


726 


1.3176 ±0.0005 


727 - 


824 


1.3158 ±0.0006 


825 - 


922 


1.3200 ±0.0006 


922 - 


1019 


1.3268 ±0.0013 



Note. — The inclination, stel- 
lar mass, and stellar radius were 
set to their best-fit values from 
Table g] for these fits. 



Table 6. Best-fit transit times 



UT Date 


N T 


T C (HJD) 


CT HJD 


(O-C) 


(O-C) 

CT H.m 


2000 April 25 


-331 


2451659.936875 


±0.000162 


0.000136 


0.84 


2000 April 28 


-330 


2451663.461599 


±0.000162 


0.000111 


0.68 


2000 May 5 


-328 


2451670.511019 


±0.000162 


0.000034 


0.21 


2000 May 12 


-326 


2451677.560259 


±0.000162 


-0.000223 


-1.37 


2003 May 3 


-18 


2452763.183042 


±0.000162 


-0.000004 


-0.02 


2003 May 31 


-10 


2452791.380849 


±0.000162 


-0.000186 


-1.15 


2003 June 25 


-3 


2452816.054642 


±0.000162 


0.000367 


2.26 


2003 July 5 





2452826.628846 


±0.000162 


0.000326 


2.01 



Note. — These are the best-fit locations for the centers of the eight eclipses 
examined in this work. The errors are based on our estimate of the size of residual 
trends in the data, which we believe to be the largest source of error. We also give 
the number of elapsed transits and O-C residuals for each eclipse. 
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Table 7. Ephemeris from this work and previous works. 



Study 


Period (days) 


T c (HJD) 


Brown et al. (2001) 


3.52480 ±0.00004 


2451659.93675 ±0.0001 


Wittenmver et al. ("2005") 


3.52474554±0.00000018 


2452854.82545 ±0.000135 


This Work 


3.52474859±0.00000038 


2452826.628521±0.000087 



